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Adiabatic dynamics of an inhomogeneous quantum phase transition: 
the case of z > 1 dynamical exponent 
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We consider an inhomogeneous quantum phase transition across a multicritical point of the XY 
quantum spin chain. This is an example of a Lifshitz transition with a dynamical exponent z = 2. 
Just like in the case 2 = 1 considered in New J. Phys. 12, 055007 (2010) when a critical front 
propagates much faster than the maximal group velocity of quasiparticles Vq, then the transition 
is effectively homogeneous: density of excitations obeys a generalized Kibble-Zurek mechanism and 
scales with the sixth root of the transition rate. However, unlike for z = 1, the inhomogeneous transi- 
tion becomes adiabatic not below Vq but a lower threshold velocity v, proportional to inhomogeneity 
of the transition, where the excitations are suppressed exponentially. Interestingly, the adiabatic 
threshold v is nonzero despite vanishing minimal group velocity of low energy quasiparticles. In the 
adiabatic regime below v the inhomogeneous transition can be used for efficient adiabatic quantum 
state preparation in a quantum simulator: the time required for the critical front to sweep across a 
chain of A'' spins adiabatically is merely linear in A'^, while the corresponding time for a homogeneous 
transition across the multicritical point scales with the sixth power of A'^. What is more, excitations 
after the adiabatic inhomogeneous transition, if any, are brushed away by the critical front to the 
end of the spin chain. 

PACS numbers: 75.10.Pq, 03.65.-w, 64.60.-i, 73.43.Nq 



I. INTRODUCTION 

A quantum phase transition is a qualitative change in 
the ground state of a quantum system when one of the 
parameters in its Hamihonian passes through a critical 
point. In a second order transition a continuous change is 
accompanied by a diverging correlation length and van- 
ishing energy gap. The vanishing gap implies that no 
matter how slowly a system is driven through the transi- 
tion its evolution cannot remain adiabatic near the crit- 
ical point. As a result, after the transition the system 
is excited to a state with a finite correlation length ^xz 
whose size shrinks with increasing rate of the transition. 
This scenario, known as Kibble-Zurek (KZ) mechanism 
(KZM), was first described in the context of finite tem- 
perature transitions [l|, |2| . Although originally motivated 
by cosmology [l|, KZM at finite temperature was con- 
firmed by numerical simulations of the time-dependent 
Ginzburg-Landau model H and successfully tested W 
experiments in liquid crystals [J] , superfluid helium 3 Q , 
both high- Tc m and low- Tc Q superconductors, and con- 
vection cells [8[. More recently, spontaneous appearance 
of vortexes during Bose-Einstein condensation driven by 
evaporative cooling was observed in Ref. [9]]. However, 
the quantum zero temperature limit, which is in many 
respects qualitatively differeiit, remained unexplored un- 
til recently, see e.g. Refs. |10l - [20| and Ref. [21] for a 
review. The recent interest is motivated in part by adia- 
batic quantum computation or adiabatic quantum state 
preparation, where one would like to cross a quantum 
critical point as adiabatically as possible, and in part by 
condensed matter physics of ultracold atoms, where it is 



easy to manipulate parameters of a Hamiltonian in time 
and which, unlike their solid state physics counterparts, 
are fairly well isolated from their environment. In fact, 
an instantaneous quench to the ferromagnetic phase in a 
spinor BEG resulted in finite-size ferromagnetic domains 
whose origin was attributed to KZM j22]. However, since 
the transition rate was formally infinite in that experi- 
ment, the KZ scaling relation between the average do- 
main size ^KZ and the quench rate has not been verified. 
The KZM argument is briefly as follows [3,[l3- When 
a transition is driven by varying a parameter g in the 
Hamiltonian across an isolated critical point at gc, then 
we can define a dimensionless distance from the critical 
point as 



9- 9c 
9c 



(1) 



When e — >■ the correlation length ^ in the ground state 
diverges as ^ '^ kl""^! and the energy gap A between 
the ground state and the first excited state vanishes as 
A ^ \e\^'^. Setting h = I from now on, a diverging A~^ ~ 
|e|^^'' is the shortest time scale on which the ground state 
can adjust adiabatically to varying e{t). A generic e(t) 
can be linearized near the critical point e = as 



e(t) 



- + 0{t% 



(2) 



where the coefficient tq is called a quench time. Assum- 
ing that the system was initially prepared in its ground 
state, its adiabatic evolution fails at an ckz when time 
tKz left to crossing the critical point equals the shortest 



time scale A -'^ on which the ground state can adjust. 
Solving this equality, we obtain 






(3) 
(4) 



From Ixz the evolution becomes impulse, i.e. the state 
does not evolve but remains frozen in the ground state at 
^KZi until —tKZ when the evolution becomes adiabatic 
again. In this way, the ground state at tKZ with a KZ 
correlation length 



e 



KZ 



^KZ 



(5) 



becomes the initial excited state for the adiabatic evolu- 
tion after —tKZ- This i^KZ determines density of quasi- 
particles excited during a phase transition in, say, one 
dimension 



e 



KZ 



(6) 



and, in general, expectation values of other operators ac- 
cording to their critical scaling dimensions. An opera- 
tor O whose expectation value scales like (O) ~ ^"^ in 
the ground state near the critical point will scale like 
(O) ~ ^KZ right after the dynamical transition. 

Note that when tq is large enough, then ^kz is small 
and the linearization in Eq. ([2]) is self-consistent: the 
KZM physics happens very close to the critical point be- 
tween —tKZ and +i.KZ- 

Recently the quantum KZ paradigm has been gener- 
alized to transitions that happen in space rather than 
in time [27H30J and to inhomogeneous transitions that 
happen in time but are inhomogeneous in space [30, l3l| . 
They are described in more detail in the following Sec- 
tions [ni and mil respectively. In an inhomogeneous tran- 
sition it is important how fast the critical front, i.e., the 
place where the control parameter e is zero moves in 
space. In Ref. [30j] we considered an example of the 
quantum Ising chain which is a representative of a uni- 
versality class with the dynamical exponent z = 1. This 
exponent means that at the critical point the dispersion 
of low energy quasiparticles is linear with a unique quasi- 
particle velocity Vq. Reference [30| demonstrates that 
when a critical front is moving with a velocity w ^ w^, 
then the transition proceeds as it were effectively homo- 
geneous and the homogeneous estimate ([6]) applies. On 
the other hand, when w <C Wq then the excitations are 
exponentially suppressed and the transition is effectively 
adiabatic. 

This picture becomes a bit more complex when z > 1 
like in the Lifshitz transition with z — 2. When z > 1 
there is a non-linear low energy quasiparticle dispersion 
at the critical point lu ^ k^ and no unique quasipar- 
ticle velocity. However, similar casuality arguments as 
for z = 1 lead to a more general conclusion that the 
transition is effectively homogeneous when v ^ Vq with 
Vq being the maximal group velocity of quasiparticles at 



the critical point (equal to the unique quasiparticle ve- 
locity when z = 1). Indeed, when the critical front is 
much faster than Vq, than even the fastest quasiparticles 
are not able to communicate any information across the 
front and the transition is effectively homogeneous. How- 
ever, we also demonstrate that the condition v <^ Vq is 
not enough for an inhomogeneous transition to become 
adiabatic when z > 1, but at the same time one does 
not need to go as far as below the minimal quasiparticle 
group velocity which is zero when z > 1. The inhomoge- 
neous transition turns out to be adiabatic below a finite 
threshold velocity v, proportional to an inhomogeneity of 
the transition, which can be obtained from a variation of 
the simple KZ argument. In this paper we both present 
the general physical arguments and support them by a 
solution of a transition across a multicritical point of the 
XY model |33| . The solution demonstrates that the gen- 
eral arguments are robust despite their simplicity, while 
the general discussion suggests that the conclusions are 
applicable beyond this particular example. 



The paper is organized as follows. In the following Sec- 
tion |lT] we provide experimental motivation and a gener- 
alized KZ argument for an inhomogeneous transition. It 
is here that we obtain the general estimate for v. In Sec- 
tion mil we do the same for a static transition in space 
and then we use the general results for a static transition 
to derive again the same estimate for v as in Section [ill 
In this way the general v is arrived at from two differ- 
ent angles. In order to test the general predictions we 
introduce the transverse field XY model in Section IIVI 
In Section |V] we review an exact solution for a homo- 
geneous transition in this model [33|, and in Subsection 
IV Al a homogenous transition in a finite chain of N spins. 
The last transition becomes adiabatic when the transi- 
tion time Tq is much longer than N^. In Section IVII 
we consider a transition in space and confirm the gen- 
eral predictions of Section iHll in the XY model. Section 
IVIII is devoted to an inhomogeneous transition in the XY 
model. In its Subsection IVII Al the effectively homoge- 
neous case oiv:^Vq is studied, and in Subsections IVII Bj 
and IVII Cj the adiabatic regime oi v <^ v. In Subsection 
IVII Bi t! is shown to be proportional to the inhomogeneity 
of the transition at the critical front. It is argued there 
that the adiabatic quantum state preparation by an in- 
homogeneous transition requires time proportional to the 
number of spins N, i.e., for large N the inhomogeneous 
transition is much faster than the adiabatic homogeneous 
transition. This conclusion is further strenghened in Sub- 
section lVII CI showing that eventual (exponentially small) 
quasiparticle excitations in the adiabatic inhomogeneous 
transition are brushed away by the critical front to the 
end of the spin chain. Finally, the concluding Section 
IVIIII provides a summary of our results. 



II. 



INHOMOGENEOUS TRANSITION 



As pointed out already in the finite temperature eon- 
text [23|, see also Ref. [24| for recent applications, in a 
realistic experiment it is difficult to make e exactly ho- 
mogeneous throughout a system. For instance, in the su- 
perfluid ^He experiments Q the transition was caused 
by neutron irradiation of helium 3. Heat released in 
each fusion event, n + ^He -> ^i?e, created a bub- 
ble of normal fluid above the superfluid critical temper- 
ature Tc. As a result of quasiparticle diffusion, the bub- 
ble was expanding and cooling with local temperature 
T{t,r) = exp{-r'^/2Dt)/{2TTDtf^^, where r is a distance 
from the center of the bubble and Z? is a diffusion coef- 
ficient. Since this T{t,r) is hottest in the center, the 
transition back to the superfiuid phase, driven by an in- 
homogeneous parameter 



e(i,r) 



T{t,r)-T, 



(7) 



proceeded from the outer to the central part of the bubble 
with a critical front rc(t), where e = 0, shrinking with a 
finite velocity v = drc/dt < 0. 

A similar scenario is generic in ultracold atom gases 
in magnetic/optical traps [23, [20| : a trapping potential 
results in inhomogeneous density of atoms p{f) and a 
critical point gc depends on atomic density p. Thus even 
a transition driven by a perfectly homogeneous g{t) is 
inhomogeneous , 



e{t,r} = 



git)-gc[pm 

9c[pir)] 



(8) 



with a surface of critical front, where e = 0, moving at a 
finite speed. 

According to the KZM, in a homogeneous symme- 
try breaking transition, a state after the transition is 
a mosaic of finite ordered domains of average size ^kz- 
Within each finite domain orientation of the order pa- 
rameter is constant but uncorrelated to orientations in 
other domains. In contrast, in an inhomogeneous sym- 
metry breaking transition (23|, the parts of the system 
that cross the critical point earlier may be able to com- 
municate their choice of orientation of the order param- 
eter to the parts that cross the transition later and bias 
them to make the same choice. Consequently, the final 
state may be correlated at a range longer than ^kz or 
even end up being the ground state, and the final den- 
sity of excited quasiparticles may be lower than the KZ 
estimate in Eq. (jS]) or even zero. 

From the point of view of testing KZM, this inhomo- 
geneous scenario, when relevant, may sound like a neg- 
ative result because an imperfect inhomogeneous transi- 
tion suppresses KZM. However, from the point of view 
of adiabatic quantum computation or adiabatic quantum 
state preparation it is the KZM itself that is a negative 
result: no matter how slow the homogeneous transition 
is there is a finite density of excitations (jH]) which decays 



only as a fractional power of transition time tq. From 
this perspective, the inhomogeneous transition may be a 
practical way to suppress KZ excitations and prepare the 
desired final ground state adiabatically. 

To estimate when the inhomogeneity may actually be 
relevant, in a similar way as in Eq. ^ and Ref. [30|, we 
linearize the parameter e(t, n) in both n and t near the 
critical front where e(i, n) = 0: 



e(t, n) « a (n — vt) 



(9) 



Here n is position in space, e.g. lattice site number, a 
is a gradient/inhomogeneity of the transition, and v is 
velocity of the critical front. When watched locally at a 
fixed n, the inhomogeneous transition in Eq. ^ appears 
to be the homogeneous transition in Eq. ^ with a local 



^Q 



1 

av 



(10) 



The part of the system where n < vt, or equivalently 
e(t, n) < 0, is already in the broken symmetry phase. An 
outcome of the transition depends on v. 

On one hand, there cannot be efficient communication 
across the critical point when it is moving faster than 
quasiparticles near the critical point; 



V ^ v„ 



(11) 



Here Vq is the maximal group velocity of quasi par ticles at 
e = or, in general, a Lieb- Robinson velocity [S^l. It is a 
constant that does not depend on the inhomogeneity a. 
In this "homogeneous regime" the inhomogeneous tran- 
sition is effectively homogeneous and the final density of 
excitations after the transition is given by Eq. ^ with 
the local tq = Xjav. 

On the other hand, KZM provides the relevant scales of 
length and time, ^kz and Ikz respectively, whose com- 
bination 1231 



e 



KZ 



^Kl 



a 1+" 



(12) 



is a relevant scale of velocity. Here we used Eqs. (|4l5p and 
Eq. ()10|) which is valid for small a. Indeed, when v <^v 
the system has enough time to adjust adiabatically on 
the relevant length ^kz and the density of excitations is 
less than in the homogeneous KZM. This is an "adiabatic 
regime" of the inhomogeneous transition. 

In general the adiabatic threshold v depends on the 
inhomogeneity a in distinction to the constant Vq. How- 
ever, in a special case of z = 1, or linear quasiparticle 
dispersion at the critical point, the two velocities are the 
same: {) ~ Wg. This (quite generic) special case was stud- 
ied recently in Refs. ;30i . i31]. In this paper we consider 
an example of a more general class of critical points with 
z > 1, or nonlinear quasiparticle dispersion, when the 
adiabatic threshold v in Eq. (|12|) scales with a positive 
power of a and, for small enough a, can be clearly dis- 
tinguished from the homogeneous threshold Vq. v <^ Vq. 



However, before we proceed with the example, in the 
next Section we rederive Eq. P^ from a shghtly different 
perspective. 



which is a ratio of the relevant length S^sp to the relevant 
time Agp. This v is the same as Eq. (IT^ . 

In the following Sections we test the general predictions 
in the XY model. 



III. 



KZM IN SPACE 



References |27H30l | considered a "phase transition 
in space" where e(n) is inhomogeneous but time- 
independent. In the same way as in Eq. ^, this pa- 
rameter can be linearized as 



e(n) sa a (n — ric) 



(13) 



near the static critical front at n = Uc where e = 0. The 
system is in the broken symmetry phase where n < ric 
and in the symmetric phase where n > ric- In the first 
"local approximation" , we would expect that the order 
parameter behaves as if the system were locally homoge- 
neous: it is nonzero for n < ric only, and when n — >■ n~ 
it tends to zero as {ric — n)^ with the critical exponent /3. 
However, this first approximation is in contradiction with 
the basic fact that the correlation (or healing) length ^ 
diverges as ^ ^ \e\^^ near the critical point and the di- 
verging ^ is the shortest length scale on which the or- 
der parameter can adjust to (or heal with) the changing 
e{n). Consequently, when approaching n^ the local ap- 
proximation {ric — n)" must break down when the local 
correlation length S, ^ [a{nc — "■)]"" equals the distance 
remaining to the critical point {ric — n). Solving this 
equality with respect to ^, we obtain 



^SP 



a i+" 



(14) 



Beginning from n — ric — —£,sp the "evolution" of the 
order parameter in n becomes "impulse" , i.e, the order 
parameter does not change until n — ric c^ +S,sp in the 
symmetric phase, where it begins to follow the local e(n) 
again and decays to zero on the same length scale of S.sp- 
A direct consequence of this "KZM in space" is that a 
non-zero order parameter penetrates into the symmetric 
phase to a depth 



Sn 



^SP 



(15) 



as if the the critical point were effectively "rounded off" 
on the length scale of £^sp- This rounding-off results also 
in a finite energy gap 



^SP 



i 



SP 



a^+" 



(16) 



in contrast to the local approximation, where we would 
expect gapless excitations near the critical point. The 
finite gap in turn should prevent excitation of the system 
even when the critical point ric in Eq. (1131) moves with a 
finite velocity: nc{t) = vt. The excitation is suppressed 
up to a threshold velocity 



e 



SP 



H^-i) 



^SP 



a i+" 



(17) 



IV. 



TRANSVERSE FIELD XY CHAIN 



The ferromagnetic transverse field XY quantum spin 
chain is 



N 



N~l 



n— 1 n— 1 

where a^^y^^ are spin-1/2 Pauli matrices. Here we 
consider a path in the parameter space {J^,jy,h) = 
(i^, ii^, 1 -I- e) parametrized by the anisotropy of the 
ferromagnetic couplings e G [—1,1]. The parameter e will 
be driven from the initial e = 1 to the final e = — 1, when 
the Hamiltonian ([T5| becomes the simple Ising chain 



N-l 



H, 



final 






(19) 



n=l 



In the thermodynamic limit a homogeneous system has 
a second order quantum phase transition (multicritical 
point) at e = 0, which separates a paramagnetic phase 
where e > from a ferromagnetic phase where e < 0. 
Non-adiabaticity of this dynamical transition will be 
quantified by average number of quasiparticle/kink ex- 
citations in the final ferromagnetic state. 

In a more general inhomogeneous system that we will 
consider in this paper it is natural to parametrize 



JX 



1 - e„ 



TV 



l + e„+i 



; f^n 



l + e« , 



where e„(t) is a continuous function of n and, in general, 
time. 

After the Jordan- Wigner transformation to spinless 
fermionic annihilation operators c„, a^ = (c„ — cjj) iSn, 
cr^ = {cn + c'l)Sn, < = 1 - 24c„, where S„ = 
n„i<„(l — 2cJj,Cm), the Hamiltonian (IT5|) becomes 



H 



E 



I - {cicn+i + h.c.) - e„+ 1 (c„+ic„ -I- h.c.) 

n 

+ (l + e„)(24c„-l) ] . (20) 

This quadratic H is diagonalized to H = J2m ^mlmlrn 
by a Bogoliubov transformation c„ = X)m=o('"n.m7"i + 
Vnmlrn) with m numerating N eigenmodes of the sta- 
tionary Bogoliubov-de Gennes equations, 

ujuf^ = 2(1 + e„)uj 

-(It e„+i )w,T+i - (1 ± e„+i)M,T-i > (21) 

with ll! > 0. Here u = w ± u. 



The Hamiltonian (|18I20I) commutes with a parity op- 
erator 



N 



N 



p = n< = n(24cn-i) 



(22) 



n— 1 n— 1 



The even parity of the initial ground state at e = 1 is 
conserved during the dynamical transition. 



V. HOMOGENEOUS TRANSITION 

When N ^ oo and e„ = e is homoge- 
neous, it is convenient to make a Fourier transform 
= ^^^^^ X]n=i ""n^**^" with pseudomomentum k £ 



(— tTjTt]. Equation ([21]) becomes 



— ( "^ ) = [a^(l + e-cosfc)-(7^esinA:l ( "^ 



with eigenfrequencies 



Wfe 



2 J (1 + e - cos fc)2 + e^ sin^ k 



(23) 



(24) 



We can expand ojq w 2|e| = |e|^'' for small e, and cjfe « 
/c^ = |fc|^ for e = and small k, to identify the critical 
exponents: 



2 , 1/ = 1/2 



(25) 



With these exponents, the general Eq. ([6]), valid for an 
isolated quantum critical point, would predict the density 
of excitations after the homogeneous transition in Eq. ([2]) 

to scale as d ~ Tq . However, as shown in Ref. |33j 
and outlined below, the multicritical nature of the critical 
point results in a different scaling ex pon ent and the KZM 
requires careful generalization here [2l|, |3J| . 

The correct scaling can be obtained by mapping our 
time-dependent problem 



i d 
2dt 



1 



cos k 



a^ — sin k 



to the standard Landau-Zener (LZ) model [3£ 



dt' 



k - 



t' 



A'ct, 



(26) 



(27) 



where t' = t + TQ{cosk — 1)/(1 + sin" fc) is shifted 
time, a — {a^ sinfc — ct^)/v 1 + sin^ k and cr^ = {a^ + 
cr^ sinA;)/v 1 + sin^ fc are two orthogonal spin compo- 
nents, r = tq /2vl -I- sin^ fc is the LZ transition time, 

and A' = 2sinfc(l — cosfc)/vl -I- sin^ fc is the minimal 
gap at the anticrossing center when i' = 0. 

In the adiabatic limit, tq ^ 1, only the long wave- 
length modes fc sa get excited with the LZ probability 



Pk 



-7rT(A') 



-TtTQk^/2 



(28) 



and the density of quasiparticle excitations after crossing 
the multicritical point is 



dk 

2^ 



Pk 



do T, 



-1/6 



(29) 



where do = 2i/6r(7/6)7r-^/6 = 0.274 and r(...) is the 
gamma function. 

The correct exponent 1/6 can be made compatible with 
the general Eq. ^ as follows. The instantaneous quasi- 
particle frequency in Eqs. (P7)) is w^ = \/{t'Y + (^')^' 
where e' = t' Jt is the relevant distance from the anti- 
crossing center. We can expand uj'y. ~ |fc|^ = |fc|^ at 
e' = and smah fc, Wq ~ |e'| = |e'|^ ^ at small e', to 
identify the exponents relevant for a homogeneous tran- 
sition as 



r' = 3 



1/3 



(30) 



Using these relevant exponents in the general Eq. © 
gives the correct exponent 1/6 in the exact Eq. (|29l) . 
Note that both z in Eq. ^, and z' in Eq. dSO]) are 
greater than 1. 

Moreover, when the transition is in space, then Eq. 
([T6)) with the exponents (|30p predicts the gap 



A 



SP 



v3/4 



(31) 



while the "canonical" exponents (1^51) give Kgp — c?!"^ . 
In Section IVII we will see that Eq. (|3ip is indeed the 
correct gap in a transition in space. 

However, depending on the choice of v or v' in Eq. (|14p 
we obtain either ^sp — oi^^^ or ^sp — ct^^^ of which only 
the former will turn out to be relevant for a transition in 
space, but even it is not a unique scale of length there. 
Nevertheless, since both z > 1 and z' > 1, we expect the 
adiabatic threshold velocity v in an inhomogeneous tran- 
sition to scale with a positive power of the inhomogeneity 
a, compare the general Eq. (IT71) . 



A. Adiabatic regime of homogeneous transition 

Since a finite chain of N spins has finite energy gap at 
the critical e = 0, the homogeneous transition becomes 
adiabatic above a finite tq when the scaling relation (P^ 
crosses over to exponential decay. 

Indeed, in a periodic chain the quasimomenta are 
quantized as k — ±-^,±-^,... to satisfy anti-periodic 
boundary conditions for the Jordan- Wigner fermions in 
the subspace of even parity. The pk = exp (^—■KTqk^/i) 
in Eq. ([^5]) is a probability to excite a pair 
of quasiparticles with momenta {k,—k). When tq 
is large enough, then only the longest wavelength 
pair (;^j— ■^) has non- negligible excitation probability 
Ptv/n — exp (— 7r^TQ/2A^^), but even this probability be- 
comes exponentially small when tq is deep enough in the 
adiabatic regime: 



TQ > 



N^ 



(32) 





FIG. 1: The critical front in Eqs. ([Ml and (gS 



The transition time required for a homogeneous transi- 
tion to become adiabatic grows with the sixth power of 
the number of spins. 

The adiabatic condition p2p can be also obtained in 
a more intuitive way from Eq. ([5|), compare similar ar- 
gument in Ref. |12| . Indeed, average distance between 

1/6 



excitations scales as ^kz ^ tJ and, consequently, the phase to a depth Sn 
tq required for the whole chain of N spins to remain 
defect- free scales as N^. This simple argument applies 
to an open chain as well. 

In Ref. [3J] a non- linear generalization e{t) = 
sign(— i) It/rgl'" of the linear quench ([2]) was proposed 
to improve adiabaticity. Indeed, a similar argument as 



X = (n-n ) a 



FIG. 2: Phase transition in space. In panel A, ferromagnetic 
magnetization X„ = {a^} is shown as a function of distance 
from the critical point n — ric for different values of the gradi- 
ent a. This ferromagnetic magnetization penetrates into the 
paramagnetic phase where n — tIc > 0. In panel B, the same 
magnetization as in panel A but rescaled by its value at the 
critical point X„^ , is shown as a function of a rescaled distance 
X = {n — nc)a^'^. The six rescaled plots collapse demonstrat- 
ing that the magnetization penetrates into the paramagnetic 



-1/3 



predicted in Eq. ([41 



a long wavelength approximation and treat n as a con- 
tinuous variable: 



in Section |T] yields dj. 



'Q 



as a generalization of 



Un±l ~ Un± -J-Un 

an 



1^ 
2dn2' 



(35) 



Eq. ©. When r > 1/zV = 1 then dr - Tq^^""' = Tq^^^ 

is much less than di ^ Tq after a linear quench for 
long enough tq. In a finite chain the non-linear quench 
becomes adiabatic when 



We also expect that the low frequency quasiparticle 
modes are localized near the critical point Uc, where we 
can use the linearization in Eq. (j34p . 
Indeed, after rescaling 



— 1/3 

n — Tir = a ' a; , u) 



a'/'n 



(36) 



TQ » N^ 
but this is still qubic in N. 



[o6) equations (PT|) become 



a 12 52 ' 

■ u 



dx^ 



+ 2x 



VI. 



TRANSITION IN SPACE IN XY MODEL 



a^/^ icjy 



d 

l + 2x— 

dx 



(37) 



In this Section we take the size of an open chain N — > 
c» to avoid boundary effects and consider a smooth static 
slant 



tanh [a[n — Uc)] ~ a{n — ric) 



(34) 



To leading order in a <C 1 their two eigenmodes of lowest 
frequency are 



r^o = 



wd ex e 2" ^ A{x) , Wq (X 



(38) 



interpolating between e — 1 and e = — 1, which is shown 
in Fig. [TJ The slant ([Ml) can be linearized near the 
critical point ric as in the general Eq. ([T3t . 

We are interested in the low frequency part of the 
quasiparticle spectrum where, presumably, we can make 



and 



r^i 



v/8r(5/4)/r(3/4) = 2.43 , 



(39) 



+ _l„l/3 2dA _ _l^l/3 2-2A{x) 

u7 GC e -^ -;-(x) , Wi ex e 2" ^ — -^—L 

dx^ ' ^ avsni 



where A{x) — Ai[2^/^a::] and Ai is the (oscillating) Airy 
function. 

The lowest energy relevant (even parity) excitation of 
the two quasiparticles 70 and 71 has a finite gap 



Agp = UIq + OJi = a ' ill 



v3/4 



(40) 



in agreement with the general Eq. (|3T|) for the choice 
z' = Xjv' = 3 in Eq. (J3 0|) 

However, the modes p8l39p do not have a unique scale 
of length. They penetrate into the paramagnetic phase, 
where n > ric and 2; > 0, to a depth 



.-1/3 



(41) 



determined by the x — ?> +cxi asymptote of the Airy func- 
tion A{x) - exp(-2y2a;3/2/3) and Eq. §^. Conse- 
quently, this Sn is also penetration depth of ferromag- 
netic magnetization into the paramagnetic phase, see Fig. 
[2j However, on the ferromagnetic side, where n < ric and 
a; < 0, the same modes (I38I39P extend to the depth 



At 



-1/2 



(42) 



limited 



by 



the Gaussian envelope e 2" ^ = 
e 2"v" ■'"^ . This envelope is damping oscillations of 
the Airy function A(x) which take place on the shorter 
scale 5n ~ a~^/^. An overall width of the modes p8l39p 
is set by the longer scale An. 

Thus in case of a finite chain size N results of this 
Section require An <C N or, cquivalently, a lower bound 



INHOMOGENEOUS TRANSITION IN XY 
MODEL 



VII. 



In this Section we make the critical front in Eq. (|34|) 
and Fig. [T] sweep from nc —t- —00 to nc —> +00 at a 
constant velocity v. 



£nit) = tanh [a{n — vt)] « a{n — vt) 



(43) 



see Fig. [TJ Near the critical point a.t n — vt the slant 
(|43)) can be linearized as in the general Eq. ([9]). 




0.01 



FIG. 3: Numerical simulations of the time-dependent 
Bogoliubov-de Gennes equations H46|) in a finite chain of 
A^ — 200 spins. Rescaled final kink density a~^'^d after 
an inhomogeneous transition in Eq. H43p is shown for dif- 
ferent values of the gradient a. In the homogeneous regime 
u 3> 2 the collapsed plots are nearly linear with a slope 0.20 
for a — 1/256 close to the predicted 1/6 in Eq. ((45)) where 



-i/«d: 



„i/6 



1/6-scaling in Eq. (|29)) to apply and a rescaled final den- 
sity of kinks to scale as 



a^i/6 d ~ z;i/6 



(45) 



with velocity of the critical front v. 

To verify this prediction we simulated time-dependent 
Bogoliubov-de Gennes equations 

*di"" = 2(l + e„)WrT 

- (lTe„+iX+i- (l±e„+iX_i • (46) 

in an open chain of N spins. Results are shown in Fig. [3] 
and in the homogeneous regime v ^ 2 they are consistent 
with the prediction (P5|). 



B. The adiabatic regime of inhomogeneous 
transition 



A. The homogeneous regime of inhomogeneous 
transition 

We can obtain quasiparticle group velocity at the crit- 
ical point e = from the dispersion (|24p . The group 
velocity is maximized for k = ±7r/2 by 

Vq = 2 . (44) 



Vq there is no causal connection across the 



When V 3> 

critical point and the inhomogeneous transition proceeds 
as if it were effectively homogeneous with a quench time 
TQ = 1/av. In this regime we expect the "homogeneous" 



In this Section we consider the adiabatic limit of small 
front velocity v when only the lowest relevant excited 
state has non-negligible excitation probability. The state 
is the even parity state occupied by the lowest two quasi- 
particles: 7o and 71. When the critical point ric is in the 
bulk of the finite lattice, then these quasiparticles are de- 
scribed by the Bogoliubov modes (I38I39P and the energy 
gap for this excitation is given by Eq. PO)) . 

In the adiabatic limit the even parity subspace of the 
Hilbert space can be truncated to an effective two-level 
system: 



m)) = a{t) |0) + b{t) |I) , 



(47) 
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FIG. 4: Parameters A and Z of the LZ problem in Eq. 
for a lattice size A*' = 400 and different values of the small 
gradient a. In panel A, the instantaneous relevant gap A = 
CJo+ciJi is shown as a function of ric. The bulk value of the gap, 
when 1 < nc < A'^, is estimated in Eq. ^^ as Agp ^ c?^^ . 
In panel B, the parameter Z in Eq. (|49|l is shown. Its bulk 



value can be estimated from Eqs. H38I39I49P as ^ ~ a 



1/4^ 

In both panels A and B, the crossover regions at both ends 
of the chain have a width An ~ a~^" set by the Gaussian 
envelope of the modes (|38l39p . The additional oscillations at 
the left end have a length scale 5n ~ q~^'^ determined by 
the oscillatory Airy function A(x) in the modes H38l39p . 



where |0) is the instantaneous ground state in the even 
subspace for an instantaneous position ric of the critical 
point and |1) = 7i7o|0) is the instantaneous first excited 
state for the same Uc- The amplitudes a, h solve a gener- 
alized Landau-Zener (LZ) problem 



dt 



ivZ 
-ivZ A 



(48) 



with initial conditions a(— oo) = 1 and fe(— oo) = 0. Here 
A = wq -l- Wi is an instantaneous gap and 



N 



Z ^ {l\- 10) = V" (Vnl.Unl 

dric ^ 

n— 1 



dric V ^no 



(49) 



Generic A(nc) and Z{nc) are shown in Figs. UA and B 
respectively. 

In the adiabatic limit of small excitation probability 
|5(cxd)P ^ 1, our numerical simulations of the LZ equa- 
tions (1^5)) are well described by a simple LZ-like formula 



H^W 



exp 



(50) 




100 

1/v 



200 



FIG. 5: In panel A, the final excitation probability |6(cx))p in 
the generalized LZ problem, defined in Eq. (|48|) and Fig. 31 as 
a function of a/v for different values of the inhomogeneity a. 
Here the solid lines are the best fits of the LZ formula (|50|) to 
the numerical data points: |&(cxd)|^ — exp (— c^). In panel B, 
the coefficients c fitted in panel A as a function of a. The solid 
line is the best fit c = 3.12 -I- 4.50a^'^ demonstrating weak 
residual dependence on a^" which becomes negligible when 
Of — >■ and c — >■ 3.12. In panel C, |6(oo)p as a function of 1/v 
for a fixed a — 2~^ and two different chain sizes A^ — 200, 400. 
In the adiabatic regime the excitation probability does not 
depend on A^ demonstrating that the excitation of the lowest 
two quasiparticles is a boundary effect. In panel D, |6(oo)p 
from the LZ model and corresponding excitation probabilities 
Po = {7o7o) and Pi = (7i7i) from the exact Bogoliubov-de 
Gennes equations (|46[) . All these three probabilities are the 
same, |&(co)| = Po = Pi, in the adiabatic regime where 
only a pair of the lowest two quasiparticles 70 and 71 can get 
excited. 



in Fig. \5\ B, there is still a weak residual dependence 
c sa 3.12-f 4.50 a^/^. However, when a < 1 then c w 3.12 
becomes independent of a as assumed in Eq. (|50p . 

Furthermore, Fig. [5] C shows that in the adiabatic 
regime the small excitation probability |6(oo)p does not 
depend on the chain size N. Not quite surprisingly, the 
excitation of the lowest two quasiparticles is a boundary 
effect determined by the behavior of A and Z in Fig. |4] 
when the critical point ric is near the ends of the chain. 



Wc can conclude that the excitation probability ^ 
is exponentially small when v <S^ a. This inequality iden 
tifics a threshold velocity 



(51) 



where c = 0(1) is a numerical pref actor. Indeed, the 
plots for different a which are collected in Fig. [S]A nearly 
collapse. The collapse is not perfect because, as shown 



when the inhomogeneous transition becomes adiabatic. 
As anticipated for 2 > 1, the adiabatic threshold t) is a 
positive power of the gradient a. 
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FIG. 6: Density distributions of kinks along a chain of A'^ = 
200 spins after an inhomogeneous transition in the adiabatic 
regime v <^ a. These distributions were obtained from nu- 
merical simulations of the exact time-dependent Bogoliubov- 
de Gennes equations H46|) with the critical front in Fig. [l] 
Residual (exponentially small) excitations are brushed away 
by the critical front to the right edge leaving behind defect- 
free bulk of the chain. The inset is a log-scale focus on the 
right edge. It shows that in the adiabatic regime all den- 
sity distributions for different v are the same up to an overall 
u-dependent amplitude set by the Landau-Zener excitation 
probability in Eq. ((50)l . 



Residual edge excitations after adiabatic 
transition 



In the adiabatic limit v <^ a only the lowest two quasi- 
particles 70 and 71 can get excited. Since the inhomo- 
geneous transition is between two gapped phases, these 
low frequency modes are localized either near the critical 
point ric when ric is in the bulk of the spin chain, or at 
one of the ends of the chain when n^ is near this end. For 
instance, the "bulk" modes in Eqs. (|38I39I) are localized 
within the distance An ~ a~^'^ from the critical point. 

Consequently, as the critical front in Fig. [1] is pass- 
ing across the chain these instantaneous low frequency 
modes follow the moving front to the right end of the 
chain. Indeed, a few generic final density distributions of 
kinks along the spin chain are shown in Fig. [51 These 
(exponentially small) probabilities of kink excitation are 
localized near the right end of the chain. Residual exci- 
tations, if any, are brushed away to the right end leaving 
behind a defect-free bulk of the chain. 



tion in an infinite chain decays with the sixth root 
of the transition time tq . 

• Consequently, the transition time required for a ho- 
mogeneous transition in a finite chain of N spins to 
be adiabatic scales like N^ . 

• Quasiparticles excited by a homogeneous transition 
are uniformly scattered along a spin chain. 

• In an inhomogeneous transition, when the criti- 
cal front propagates much faster than the maximal 
quasiparticle group velocity w^ = 2, then the tran- 
sition is effectively homogeneous. 

• When a critical front moves much slower than the 
adiabatic threshold velocity w ~ a, then average 
number of excitations is exponentially suppressed 
and does not depend on N. 

• Quasiparticles excited in the adiabatic regime, if 
any, are brushed away by the critical front to the 
end of a spin chain leaving behind defect-free bulk 
of the chain. 

• The minimal time required for the adiabatic inho- 
mogeneous transition to be completed is N/v. It 
is a mere linear function of N instead of the N^ in 
the homogeneous case. 

The main difference between the case 2; = 1 in Ref. fSO] 
and z > 1 considered here is that in the former case the 
adiabatic threshold velocity v is the same as the homo- 
geneous threshold velocity Vq, v = Vq when z = 1, while 
in the latter case the adiabatic threshold is less than the 
homogeneous one, v < Vq when z > 1, and it is pro- 
portional to the inhomogeneity of the transition. When 
z > 1 the KZ adiabatic threshold v is nonzero despite 
vanishing minimal quasiparticle group velocity. 

Putting the results for z > 1 together with those of 
Ref. [30| for z = 1, we can conclude that for large N 
an inhomogeneous transition is a more efficient method 
of adiabatic quantum state preparation than a straight- 
forward homogeneous transition. Not only the time re- 
quired for an adiabatic transition is much shorter, but 
also any residual (exponentially small) excitations are 
brushed away to the end of the spin chain. 



VIII. 



CONCLUSION 
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